Is there ERE in bright nebulae? 
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' Abstract. In 1986 Witt & Schild studied the optical emission of selected bright nebulae. 
Comparison of the data to a model led the authors to conclude the nebulae have too high a 
OO ' level of /-band emission to be explained solely by the scattering of starlight. The excess is 
explained as a broad band / emission which originates in the nebulae. This phcnomcnum 
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is known as Extended Red Emission (ERE). 

Given Witt & Schild's data, the model they employed, and the estimate of obser- 
vational errors, I've reached a different conclusion: Witt & Schild's observations are 
compatible with classical scattering and therefore cannot be used to prove the existence 
of ERE. 
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1. Introduction 



In 1986, Witt and Schild (198(;, WS hereafter) presented CCD photometric observations 
of 15 nebulae, amongst which are the brightest nebulae of the Galaxy, all illuminated by 
close B stars. 

The data set assembled by WS is of great interest since it fosters better understanding 
of the properties of interstellar grains and of the structure of the nebulae. The dependence 
of the surface brightness of a nebula with distance to the star, the estimate of the surface 
brightness per unit radiation field, and the color of the nebulae, narrow the search for 
better knowledge of grain properties and will indicate whether or not scattering alone 
accounts for all nebulae emission. 

The surface brightness at wavelength A, S\, of a nebula assumed to scatter the light 
of a nearby star is proportional to the starlight flux at the cloud location. Wc have 
no direct way of knowing but the ratio / Sx^ of two optical bands is proportional 
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to / F^^ , being the flux of the star measured on earth and corrected for reddening. 
S\ is an observational value and F^ can be estimated from the observed flux of the 
illuminating star and the extinction Ay in the direction of the star. 

AC"(Ai, A2) = \og{S\-^ / S\r, X F^^/F^_^) can be known with sufflcient accuracy to be 
compared to simple models, such as the ones proposed by Witt (1985) or by Zagury, 
Boulanger and Banchet ( 1999| ). AC°(Ai, A2) depends exclusively on the optical depth of 
the scattering medium and on the properties of the grains. It may also depend on the 
structure of the medium if regions of different optical depths are mixed in the beam of 
the observation. 

The WS paper compares the colors AC°(i?, V) and AC''(y,/) of the nebulae to the 
model proposed in Witt ( 1985| ). The conclusion is summarized by flgure 14 of WS' paper, 
which is reproduced here (figure Q). WS attribute the difference between the model and 
the observations to a non-scattering emission process in the / band. This Extended Red 
Emission (ERE) is presumed to be luminescence from a particular type of interstellar 
grain. 

Before agreeing with WS' conclusion different questions must be clarifled. 

Why do most of the observational points of the WS figure 14 follow a curve which is 
parallel to the model? 

How does the reddening of the stars affect AC" and modify the relative positions of 
the observation values and the model on the color-color plot? 



According to Witt ( 1985 ), no model can be perfect, and models of reflection nebulae 
usually meet with only limited success. The Witt model, which represents nebulae as 
homogenous gas masses of defined geometry, may incorrectly represent the nebulae. Will 
small scale structure, for instance, modify the color forecast of an observation from the 
one predicted for a homogenous nebula? Are the limits of the model consistent with the 
high column densities which likely prevail in the nebulae? 

How do observational and data reduction errors influence WS' conclusion? Variations 
of surface brightness with angular distance 9 to the star, which roughly follows a 9^^ will 
be used to prove the abnormal behaviors of the I surface brightness of some nebulae. A 
major, and extremely difficult to correct, source of error is the gradient due to starlight 
diffusion in the earth's atmosphere which usually accompanies observations of nebulae 



at close distances to a bright star (see Zagury, Boulanger and Banchet, 1999). 

Attempts to answer these questions will be found in the following sections. In section^, 
WS data and argumentation are reviewed. In section |3| AC'^{B,V) and AC°(V,/) are 
extracted from the data and compared to the Witt model. The influences of column 
densities higher than permitted by the low column density approximation of the Witt 
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Fig. 1. left: WS figure 15, AC{B,V) versus AC{V,I). The line drawn corresponds to 
AC'^{B,V) versus AC°(V,/) for a homogenous mass of gas of varying optical depth, 
computed in the approximation of single scattering, right: observed AC'^{B,V) versus 
AC°{V, I). Calculated from E{B ~ V) given in WS. This plot does not figure in WS. 

model and the structure of the nebulae will be discussed to reach an understanding of 
the differences between the model and the observations. 

The general idea which will emerge from this study is that scattering can explain all 
the nebulae emission. It renders unnecessary the introduction of ERE carriers to explain 
the optical emission of this sample of nebulae. 

2. Data and WS interpretation 

For 3 to 25 chosen positions in each field, WS give the distance to the star, dg, and 
\og{S\/ F^), where S\ and are the nebula surface brightness and the star flux (not 
corrected for reddening), measured on earth, in the V, B, R, I bands. Because of the 
proximity of the star (positions are between 20" to 150" to the star), WS had to subtract 
a gradient due to starlight scattering in the earth's atmosphere. 

I used WS data to reproduce their figure 15 (figure |l| left). The plot is AC{V, I) versus 
AC{B, V) with AC(Ai, A2) = log(5A,/FXj - \og{S^jFl). 

The plain line in the plot is the variation with optical depth of AC{V, I) and 
AC{B, V), and assumes no reddening of the illuminating star. It was calculated from the 
model described in Zagury et al. ( |l999|) whose application is restricted to single scattering 
and is in agreement with the Witt model. Both models assume the star and the nebula 
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to be equally reddened by forward material -if any. Witt's model (see the equation 2 



of Witt 1985) approximates the scattered part of the extincted light by g((p)w(l — e^"^) 
where uj is the albedo of dust grains, g the phase function, ip the angle of scattering, and 
T the optical depth of the medium. The model applies to small optical depth mediums 
where single scattering dominates. 

Most of the observed points are below the model line, which WS interpret as / values 
too large to be the result of scattering. They conclude that ERE accounts for 30% to 
50% of the / band surface brightness. 

3. Discussion 

3.1. Correction for the extinction of the stars 
I will introduce: 

AC"(Ai,A2) =l0g(5Ai/FAj -l0g(^Aji^Aj (1) 

Fx is the star flux at wavelength A, measured on earth and corrected for reddening. 

AC° so defined is a possible definition for the color of a nebula. A significant difference 
exists between the traditional meaning of color and AC*': when increasing the optical 
depth of a low density medium, B — V, B — I and V ~ I colors increase while AC'^{B, V), 
AC"(B,/), AC°{V,I) decrease. 

AC"(y, /) and AC^{B,V) anticipated values for a medium with little reddening 
are AC"{V,I)max = 0.34 and AC'^{B,V).,nax = 0.13 {Ab/Av = 1.34 and Aj/Av = 



0.45, Cardelli et al. 1989 ). Extinction effects in mediums with increasing optical depths 
can only decrease (redden) AC°(y,/) and AC'^{B,V), so that one should always have: 
AC"(y,/) < AC'^{V,I)max, and AC"(B,y) < AC'^{B,V)max, if scattering alone is 
considered. 

The model curve of figure |l| is the theoretical color-color plot (AC°(F,/), AC^{B,V)) 
for a medium of low optical depth. The model assumes that starlight is not reddened 
between the star and the nebula. Starlight extinction between the star and the scattering 
volume will redden the visible emission of the nebula and move the model curve down 
and to the left, thus approaching the observations. 

AC°{V, I) and AC"{B, V) can be estimated from AC{V, I) and AC{B, I) in the fol- 
lowing manner. At optical wavelengths. Ax is a linear function of 1/A (Cardelli, Clayton & 
Mathis |1989| , Rieke & Lebofsky [l985| ). Hence: 



1_ \ ^ 

A R A V ^ 



2.2E{B^V)(^-^) (2) 
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Use of relation || in the definition of AC° gives: 

AC"(Ai,A2) - AC(Ai,A2)-0.9i?(B~T/)(^^-^) (3) 

AC°{B,V) ^ AC{B,V) ~OAE{B ~V) (4) 

AC°{V,I) = AC(V,I) -0.6E{B -V) (5) 

Figure |l|, left, plots the observed values of AC{V,I) and AC{B,I) for the nebulae WS 
have considered. Most points follow a curve parallel to the model curve but shifted down 
from it. 

A few points, mainly observations of NGC2068, have AC{V,I) greater than 
AC'^ {V, I)max- Roughly half of AC{B,V) values of the plot are significantly above 
0.13 mag. This is to be attributed to the reddening of the stars which decreases 
AC^{B, V) and AC°(V, /) according to equations ^ and ||. Correction for the reddening 
of the stars will shift all the observations down by 0.6E{B ~ V) mag and to the left by 
OAE{B - V) mag. E{B - V) is given in table 1 of WS. 

The color-color diagram (AC°(y, /), AC°(y, /)) is represented on the right hand plot 
of figure |l|. Most nebulae observations are now scattered close to the model curve. 

Different reasons may explain departures from the model line and the tendency for 
the nebulae to remain under the curve. Error bars explain part of the difference between 
the model and the observations. The Witt model is restricted to low column densities 
since it requires single scattering and low optical depth. What is the effect of an in- 
crease of column density? The nebulae of the Witt model must be homogenous enough 
for the medium sampled by the beam of the observation to be described by one value 
of the optical depth {tq in Witt's model). How will small scale structure modify WS' 
conclusions? 

3. 2. The limits of the Witt model 

The limit of validity of the model curve (figure |l|) is that single scattering dominates 
for each of the optical bands. It is satisfied if the scattering volumes are proven to be 
homogenous (can be described by a single optical depth) and have a low Ay value. 



According to Witt ( 1985 ) the model is valid for Ay ^ t < Tum — 0.6. Since this 
condition must be satisfied for all the optical bands, it implies: Ay ry < 0.6, hence 
Ai - O.bAv < 0.3. 

WS' representation of the medium surrounding the stars (§ III b in WS), by a ho- 
mogenous spherical nebula of radial optical depth Ay ~ 0.5 mag, is within the limits of 
the single scattering model but is questionable. The selection criteria of nebulae with high 
surface brightnesses in largely extincted regions, some of which are starforming regions. 
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supports the idea of high cohimn densities and perhaps a wide range of structures with 
different optical depths along each line of sight. 

High column densities in the nebulae are also indicated by the extinction in the 
direction of the stars. Ay in the stars direction should be the radial optical depth to use 
in the Witt model. Nearly all the stars in the sample have large reddening. E{B — V) for 
all stars except one, the illuminating star of IC426, is greater than 0.2, with a mean value 
- 0.7 (WS, table 1). With a minimum Ry = Av/E{B -V) oi'i (Cardelh et al. |l98|). 



the Av = 1.08r* value found for the nebulae will range from 0.6 to 4.2, with an average 
value of 2, far above and above the radial optical depths adopted in WS. 

Observations do not agree with WS' concept of the nebulae. Column densities must 
be higher than Ay — 0.5. The effect on the color-color plot of such high column densities, 



out of the range of validity of the single scattering model, is investigated in section 3.3. 

A clumpy medium made of cells of different column densities, and much smaller 
than the resolution of the observations, may be an alternative representation of the 



homogenous mass of gas adopted in WS. It will modify (section 3.4) the color of the 



nebulae and can explain the high surface brightnesses observed in all the optical bands. 
3.3. Effect of an increase of the column density on the model curve 

If the column density of the medium which is observed is increased, the single scattering 
approximation will stop being valid in the B and in the V bands before the / band. At 
low column density [Ay ^ 1), scattering in the / band will be 3 times less (Cardelli 



et al. 1989) than in the B band. Increasing the column density will sec the domination of 
absorption effects in the B band while scattering in the / band will become important. 
For Ay ^ 3, we have Ai ^ 1.35 and As ~ 4. Given these values, scattering will be 
minimal in the V and B bands, most of the B and V light will be absorbed. It is close 
to its maximum in the / band. B and V surface brightnesses will be close to and / 
surface brightness close to its maximum. 

Hence, the net effect of an increase of column density will be to shift the model points 
of figure down, as observed. The points are also shifted to the left, but the effect is less 
important (section 3T, equations ^and||). 



The model curve of figure |^ is a limit which separates two regions. All observations 
should lie in the region under the curve. The few points above the model curve will be 
explained by the error margin of the observations. 



3.4^. Effect of the small scale structure 

The ambiguity between the high column densities which exist in the nebulae and the 
observed high surface brightnesses in all bands can be removed if the emission in the 
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Fig. 2. Plot of log(S'A/Fx) versus distance to the illuminating star for the B, V, R, 
I bands and for 3 nebulae. The / band values of IC435 (left plot, asterix) are nearly 
constant and do not follow the variations of the other band, as in most other fields. 
A plausible interpretation is the gradient subtracted from 5a, which may introduce a 
variable offset. 

various bands is not exactly produced by the same interstellar regions on the same line 
of sight. High / band surface brightness arises from regions with relatively high optical 
depth {Ay > 1) which will absorb V and B starlight and produce little emission in those 
bands. The B and V band emissions should come from regions of the same line of sight, 
with lower Ay, as in the WS model. Those regions will give little scattering in the / 
band. If this is the case, one cannot expect the observations to match the uniform mass 
of gas model. 

The contribution of structures with high Ay to the emission in the /-band will con- 
tribute to the drop in AC'^{V, I) observed for some points in the right plot of figure 0. 

The probability that the observed nebulae are structured at small scales is strength- 
ened by the beam size of WS observations. WS images' resolution is not better than 
10" and the average distance to the nebulae is, according to WS' table I, 600 ± 300 pc. 
Thus, the spatial resolution of the observations is at least 0.03 pc, a large value (3 orders 
of magnitude) compared to the possible interstellar clouds scale (~ a few tens of A.U., 
Falgarone et al. |199S|, Zagury et al. |1999D. 



3.5. Error bars 

Most observations (figure |l|, right) lie under the model line and at close distance to it, 
which can be justified either by higher column densities than permitted by the model or 
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by the small scale structure of the medium. Error bars will explain the extreme position 
on the plot of some observations. 

The uncertainty of E{B — V) introduces an error on AC°{B,V) and on AC^{V,I). 
The error on E{B ~ V) is probably less than 0.1 and may produce a similar error in 
AC"{B,V) and AC"{V,I). 

Error bars on the observed values, AC{B, V) and AC{V,I), may be larger than WS 
estimates due to the atmospheric gradient substracted from each observation. The same 



problem was dealt with in Zagury et al. ( 199£ ): a precise subtraction of the atmospheric 
gradient was found to be impossible to arrive at. 



3.5.1. NGC1333 and NGC2247 

The singular positions above the model line of NGC1333 and of some of NGC2247 ob- 
servations must be due to an error in the data. For NGC2247, the error, given in WS, 
- ±0.04 to ± 0.12 for AC{B, V) and AC{V, I), is enough to bring all points back under 
the model line. 

Errors given for NGC1333 are lower, of order ±0.04, and may have been slightly under- 
estimated. Along with a possible overestimation of the E{B~V) value of BD+30°549, the 
illuminating star of NGC1333, error bars also explain the extreme positions of NGC1333 
observations in the right plot of figure n. 



3.5.2. IC435, NGC2245 and NGC6914a 

Three nebulae, IC435, NGC2245, NGC6914a are far under the model curve. Some of 
these points have AC'^{B,V) > AC''\B,V)max which must indicate error bars larger 
than WS' estimates. 

Figure ||displays log(S'A/-F\) versus logdg for some stars. For most nebulae the surface 
brightness decreases when moving away from the star, in all the optical bands. The de- 
crease follows a l/d^g law for the fields NGC1788, IC446, NGC2247, NGC2327, NGC7023, 
and Ced201 (figure ^, right plot). Similar variations of log(S'A/F\) versus logde are ob- 
served in most fields and all optical bands. Important differences also arise, especially in 
the / band. 

A particularly neat example is IC435 where B, V, R surface brightnesses decrease as 
1/dg, while / surface brightness remains constant (figure ||, left plot). hog{S\/FJ^) values 
for IC345 are nearly equal in the B, V, R bands (figure |l|) at all positions. The peculiar 
behaviour of IC435 in the /-band (figure |^, left) is most probably explained by a large 
error in / surface brightness. The error on AC{B, V) and AC(y,/) for IC435, estimated 
by WS, is of order ±0.1. An error of ±0.2 to ±0.3 seems to be more likely for AC{V, I). 
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Like IC435, the 7-band variations (not represented here) of NGC2245 and NGC6914a 
do not follow the B, V, and R band variations. Errors for these two stars, given in WS, 
are substantially higher than for IC435. They reach 0.2, a value comparable to the one 
infered for IC435. 

4. Conclusion 

In this paper I tried to show that WS analysis of the emission of bright nebulae is not 
conclusive. The observations are compatible with starlight scattering as the only process 
involved. 

The first part of WS' paper is dedicated to the comparison of the color of the nebulae 
with Witt's model, which applies to low column density nebulae where single scattering 
in all optical bands dominates. The curve which represents this model in a (AC°(B, V), 
AC°(V,7)) plot is to be considered as a limit under which all observations should be. 

Correction for the reddening of the stars will scatter most observations close to the 
model line, where they have a tendency to remain under the line, as expected. 

Three reasons explain the remaining differences between the model and the observa- 
tions. 

A deeper analysis of the error margin of the observations explains the largest differ- 
ences between the Witt model and WS' observations. The main source of error is the 
gradient of starlight scattered in the earth's atmosphere which needs to be subtracted 
from the observations. 

The general drop of the observations under the model curve can be attributed to 
column densities higher than tolerated by the single scattering model. Fitting the obser- 
vations to the model in the B or the V band, as is done in WS, will give low Aj and 
low I surface brightness values for the nebulae, incongruous with the large reddening of 
the stars. The Ay values in the direction of the stars are between 1 and 4, with a mean 
value of 2, above the maximum optical depth (~ 0.6) supported by the model. 

High surface brightness in the I-band does not prove ERE, but indicates optical depths 
in the B and V bands higher than can be supported by the model. The expected position 
of these observations on the color-color plot will be under the curve, as it is observed. 

The small scale structure of the nebulae will also heighten differences between the 
Witt model and the observations. A structured medium permits high surface brightnesses 
in all bands since the bulk of the emission at two different wavelengths need not arise from 
the same interstellar structures. Small scale structure in the sample of nebulae presented 
in WS is probable because of the large beam of the observations and the high column 
densities in the nebulae. 



10 Frederic Zagury: ERE in bright nebulae 

The Witt model, which represents the nebulae as spheres of constant density centered 
on the star, never fits the V band observations variation with distance to the central star. 
But up to what point is it reasonable to expect the observations to match the model? 
When the observations do not match the representation of a nebula proposed in WS 
shouldn't the model be questioned first, before concluding that interstellar grains have 
special properties which vary from cloud to cloud? 

The many different aspects discussed in this paper indicate that WS' observations can 
be explained by scattering effects only. While it certainly is possible to interpret the data 
with excess of red emission, this interpretation requires a constrained and limited model 
which is probably a poor representation of reality and whose reliability the authors have 
yet to prove. 
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